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Abstract. 

We report on a photometric study of a sample of 22 disk galaxies in the Hubble Deep Field South NICMOS parallel 
field. The redshift range of the galaxies is z = 0.5 — 2.6. We use deep NICMOS J and H band and STIS open 
mode images, taken as part of the HDF-S project, to construct rest-frame B-profiles and (U — V) color profiles 
of the galaxies. Before fitting isophotes, images are deconvolved with PSF. Derived surface brightness profiles are 
approximated by Sersic luminosity distribution. Significantly large population of disks cannot be represented by 
an exponential disk, but this can be well done by Sersic law, if n < 1. This might be the same phenomenon which 
has earlier been referred to as truncation of disks. Parameter n does not vary significantly with redshift. Galactic 
sizes decrease with redshift as r e (z)/r e (0) = 1 — 0.26 2 . The rest frame (U — V) color shows a clear decrease at 
2 ~ 2, concordantly with the understanding of more intense star formation at earlier epochs. Color gradients 
A(U — V)/Ar are small and roughly constant at 2 < 2. At 2 > 2, dominantly positive gradients appear, possibly 
indicating centrally concentrated star-formation. On the basis of (U — V) color and chemical evolution models, 
the disks observed at 2 ~ 2.5 have formed between 2 = 3.5 — 7. Scale radii r e of the galaxies correlate with the 
scale surface brightnesses for the sample. None of the studied parameters shows clear dependence on absolute 
B luminosity for the galaxies. 

Key words, galaxies: photometry - galaxies: fundamental parameters - galaxies: high redshift - galaxies: spiral - 
galaxies: structure 


1. Introduction 

Comparison of the structure of galactic disks at inter¬ 
mediate and high redshifts with the disks resulting from 
N-body simulations allows us to test the corresponding 
algorithms and to specify physical processes involved in 
disk formation. General properties of both elliptical and 
disk galaxies are rather well reproduced by modeling (e.g. 
White & Rees llf)7%l Navarro & White IT^Ml Abadi et al. 
120031 Bell et al. 120031 Governato et al. 120041 Robertson et 
al. 12223 Nagamine et al. 12005tl . The initial shortcomings 
have been corrected by now and it has become possible to 
model the formation of extended, rotationally supported 
disks. On the other hand, good models for calculating the 
formation and evolution of galaxies involve a great number 
of free parameters (see e.g. Cole et al. ESP. 

To test and constrain the simulations carried out at 
increasing resolutions, detailed observational studies of 
individual high redshift galaxies are indispensable. In a 
number of papers, structural parameters of galaxies have 
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been statistically studied at a wide range of redshifts. 
Luminosity function has been derived out to 2 ~ 1.5 —2 by 
Chen et al. cnna and Ilbert et al. C2E3; to 2 ~ 3.5 by 
Poli et al. TO and Giallongo et al. GDQ3; to 2 ~ 5 by 
Gabasch et al. rcnna . General morphologic studies and 
statistics out to redshift 2 « 3 and even further have 
been conducted by Bunker et al. in™ . van den Bergh 
et al. 1200011 . Corbin et al. I|2000l) . Bouwens et al. 1200411 . 
Cassata et al. (EH3 and many others. 

Due to hierarchical clustering, galactic sizes should 
grow with time (e.g. Fall & Efstathiou EMU- An anal¬ 
ysis of galactic size distribution allows us to constrain 
certain parameters in galaxy evolution simulations (feed¬ 
back etc., see Cole et al. Enn@. The observations do not 
clearly confirm this prediction as yet. For example, galaxy 
size evolution with redshift has been studied by Ferguson 
et al. 1200411 and Trujillo et al. ( 20041 BUTTS ). They show 
that, at a given luminosity, decrease of the sizes of galax¬ 
ies with redshift is clearly present. On the other hand, 
Ravindranath et al. ITOl and Cassata et al. H2005I) do 
not detect significant evolution. Serious sources of contro- 
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versies are selection effects and difficulties in conducting 
photometric analysis of galaxies at high redshifts, both 
caused mainly by cosmologic dimming of surface bright¬ 
ness by a factor of (1 + z) 4 . 

While there is a growing inflow of data concerning gen¬ 
eral statistical properties of high redshift galaxies, just a 
few papers have concentrated on the internal properties 
of distant disk galaxies, e.g. kinematics, detailed surface 
brightness and color distribution. For example, rotation 
curves for distant disk galaxies (z < 1) have been mea¬ 
sured by Vogt et al. ( 119561, 1199711 . Erb et al. (120031) and 
Bohm et al. (1200411 . Moth & Elston (1200211 have studied 
HDF-N field and constructed rest-frame (UV As — U 300 ) 
color profiles and rest-frame B surface brightness profiles 
for 83 galaxies at z = 0.5 — 3.5. They report color gradient 
A(t/V 218 — f7 3 oo)/A?' rising with redshift, which suggests 
that star formation has shifted inwards during the evolu¬ 
tion of galaxies. 

In Tamm & Tenjes 1)20031 200 5). we combined lumi¬ 
nosity distribution with rotation curves to construct self- 
consistent mass distribution models for disk galaxies out 
to z « 1. We noted that for all the 7 galaxies studied, 
the usual exponential disk assumption gave poor fit to 
the luminosity profiles; a more steep cut-off at outer radii, 
provided by the Sersic CMEi index n < 1, was needed. 
Was this because of selection effects, evolution of disk pa¬ 
rameters with time or anything else? With this question 
in mind, in the present paper we concentrate on the lu¬ 
minosity and color profiles of disk galaxies at a redshift 
range as wide as possible, in an attempt to detect evolu¬ 
tionary effects of several disk parameters with a particular 
emphasis on the shape of the profiles. To achieve this, our 
main task is to acquire surface brightness profiles as well 
as color profiles for a sample of disk galaxies, introducing 
as few selection effects as possible. 

We have chosen a small field in the southern sky, the 
Hubble Deep Field (HDF) South NICMOS parallel field 
for the study. In addition to the information offered by 
the multi-color observations of this field, the choice of the 
southern sky provides comparison with the well-studied 
HDF-N region. 

In the present work we take Ho = 65 km s _1 Mpc _1 
and Ho = 1, fl m = 0.3, Ha = 0.7. 

2. Observations and sample selection 

To avoid generation of artificial evolutionary trends, 
galaxies at different redshifts should be observed in the 
same rest-frame waveband. This precaution eliminates 
possible effects of “morphological k-correction”, i.e. the 
dependence of the photometric properties of a galaxy on 
the rest-frame waveband. This phenomenon is especially 
frequent in the case of mid-type galaxies (Papovich et al. 
I2M31) . For comparison with near-by galaxies, rest frame 
optical passband would be the best choice. To study rest- 
frame optical properties of galaxies at high redshifts, ob¬ 
servations made in near infrared are required. Presently, 
the best available resolution in near infrared is offered by 


the NICMOS camera aboard the Hubble Space Telescope 
(HST). However, the imaging properties of the NICMOS 
camera are still far behind those of the HST optical cam¬ 
eras. 

During HDF South observations (Williams et a.l. 1200011 . 
the NICMOS camera was pointing in a slightly different 
direction, and in this way, a parallel field was created with 
NIC-3 camera through the J, H and K broad-band in¬ 
frared filters. This field is also covered by WFPC2 cam¬ 
era / filter observations of Flanking Field 9 (Lucas et al. 
12003 ). Additionally, the NICMOS parallel field has been 
observed with the HST STIS camera using the open mode, 
i.e. without any filter, which gives a very broad bandwidth 
with the central wavelength matching that of the stan¬ 
dard U-band. Regrettably, the / and K band observations 
could not be used here - the former suffers from insuffi¬ 
cient depth and the latter has too high noise and back¬ 
ground level (the K images were taken during ’’bright” 
time, i.e. when the telescope was pointed near the bright 
limb of the earth). 

Observations conducted in the three remaining pass- 
bands - V, J and H - have been used in the present study. 
The basic properties of these observations and the final 
images are given in Tabled 

We have used the fully calibrated, combined and 
dithered images of the HDF-S NICMOS parallel field, 
available via the homepage of the Space Telescope Science 
Institute. The lower limit of detectable surface brightness 
on these exposures is approximately 27.5 mag/arcsec 2 . 

Table □ shows that on NICMOS exposures the full- 
weight half-maximum (FWHM) of the point-spread func¬ 
tion (PSF) is 0.23", causing serious distortion of the galac¬ 
tic images. The width of the PSF as well as its slightly 
non-circular shape cause uncertainties in deconvolution of 
the images. 

Unfortunately, spectroscopic redshifts of the galaxies 
in the HDF-S NICMOS parallel field have not been mea¬ 
sured. Photometric redshifts, based on 9-band measure¬ 
ments, have been calculated by Yahata et al. rcnnra . 

Galaxies with redshift z > 0.5 were selected for the 
further study. Discrimination between Hubble types is not 
always a simple task at high redshifts. Many objects may 
be in transition stages (protodisks or protospheroids) even 
at 1 < z < 2 (Conselice et al. 120041 ). Morphological struc¬ 
ture of these objects was inspected visually on the ba¬ 
sis of H -band exposures, all galaxies suspected for be¬ 
ing disks were included in the sample. In some cases, also 
the high-resolution STIS U-band images were checked (see 
Table 1 for the comparison of PSFs and projected pixel 
sizes). Final confirmation for the sample to consist of disk 
galaxies comes from the luminosity profiles, which all ex¬ 
hibit Sersic index n < 2 and should thus have a late-type 
morphology (Andredakis et al. Il9931 Ravindranath et al. 

TO . 

Galaxies showing significant asymmetry or irregular 
shape were excluded, because fitting ellipses to their 
isophotes would give rather scattered light profiles (light 
distribution is very sensitive to the galactic center posi- 
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Fig. 1. Galaxies studied in the present paper as they appear on STIS open mode (~ V) and NICMOS J passband 
images of the HDF-S NICMOS field. For clarity, the studied galaxy has been encircled on some images. The angular 
size of each image is 4.5 x 4.5 arcseconds. x and y coordinates of NICMOS images are used as galactic names; z is 
photometric redshift. 
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Fig. 2. Rest-frame R-band luminosity profiles with Sersic fits, and rest-frame (U — V ) color profiles with linear fits. 
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Fig. 2. (Continued) 
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Table 1 . General information about the HDF-South NICMOS field observations 


Camera & filter 

Central 

wavelength (nm) 

Bandwidth 
(FWHM, nm) 

Total exposure 
time (s) 

PSF width 
(FWHM, ") 

Projected pixel 
angular size (") 

STIS open mode 

550 

441 

25900 

0.065 

0.025 

NICMOS FI 10W 

1132 

588 

108539 

0.23 

0.075 

NICMOS F160W 

1608 

399 

128441 

0.23 

0.075 


tion) and their interpretation would not be straightfor¬ 
ward in the context of the present models and subsequent 
analysis. Thus possible starburst and interacting galaxies 
were rejected as photometrically incomparable to regular 
disks. 

These selection criteria finally set a redshift limit at 
z = 2.6, beyond which no disk galaxies could be distin¬ 
guished with acceptable confidence. A final sample of 22 
galaxies qualified for further analysis. 

The choice of passbands described above allows us 
to determine rest-frame optical luminosity profiles, using 
STIS observations for redshifts z < 1.0 and NICMOS J 
and H observations for 1.1 < z < 2.0 and z > 2.1, respec¬ 
tively; the mean central rest-frame wavelength thereby be¬ 
comes 420 nm, corresponding to Johnson B filter. Color 
information can be obtained, using STIS and NICMOS J 
observations at z < 1.1, and NICMOS J and H bands at 
z > 1A (no galaxies were found at 1.1 < z < 1.4). The 
mean central rest-frame wavelengths thus become 350 nm 
and 580 nm, allowing the derivation of (U — V) color dis¬ 
tribution. Note that here STIS observations are used for 
U and B rest-waveband photometry for the same galaxies. 
This can be justified by the very wide “passband” of STIS 
open mode, depending only on the detector sensitivity. A 
drawback of using such a wide wavelength range lies in the 
danger of suppressing possible color-features and trends, 
but it does not introduce or artificially amplify the evolu¬ 
tionary effects. 

The STIS and NICMOS J-band images of the galax¬ 
ies are presented in Figni STIS images are designated as 
V; x and y pixel coordinates of the galaxies on NICMOS 
images are used as galactic names; z is photometric red- 
shift. Where necessary, the galaxies of the present sample 
have been encircled to avoid confusion. Some of the sam¬ 
ple galaxies lie at the edge of the field of view of the STIS 
camera and the galaxy 380-1027 falls just outside it, thus 
no color-information could be acquired in the latter case. 

3. Data reduction and photometry 

To determine absolute magnitudes and especially to gen¬ 
erate color profiles, it is vital to precisely estimate the 
background level for each galaxy. We measured back¬ 
ground level at 2 — 4 empty-looking fields around each 
galaxy and found the uniformly subtracted background 
levels of both STIS and NICMOS final images to be 
slightly over-estimated. Individual background estimates 
for each galaxy were necessary due to considerable varia¬ 


tions within each passband image, ranging from —0.00015 
to 0 counts per second in the case of STIS and from 
—0.00005 to 0.00002 counts per second in the case of 
NICMOS images. 

Before fitting ellipses to the galactic images with the 
STSDAS task ELLIPSE, we deconvolved images with 
PSF. Correct deconvolution is very important, if one 
wishes to conduct analysis of the luminosity and color pro¬ 
files of faint objects. The bulge component of a disk galaxy 
often acts as a point source at high redshifts; wrong esti¬ 
mating of the spreading of its flux may cause significant 
drift of the parameters of the disk component. 

In the case of the NICMOS images, the model PSFs 
were created using stellar images found in the same field. 
In the case of STIS, no unsaturated stellar images were 
found and the model PSF was created using the TinyTim 
package. 

Ellipses were interactively fitted to the isophotes. Any 
object, which did not seem to be part of a given galaxy, 
was masked. Both fixed and free wandering values were 
tested for ellipse centers; the resulting differences in sur¬ 
face brightness profiles are included in the uncertainties of 
the profiles in Fig. [3 Ellipticity and position angle were 
kept fixed according to the outer isophotes of the galaxies. 

For rest-frame photometry we transformed the ob¬ 
served passbands into rest-frame standard passbands, 
which needed minimal k-correction after redshifting. B- 
profiles and (U — V ) color profiles were found to be most 
suitable (see Sect. 2). 

No reliable k-corrections existed for the transforma¬ 
tions from neither the STIS open mode nor the NICMOS 
passbands and we had to calculate them. In our calcula¬ 
tions we relied on the synthetic spectra of redshifted Sb 
galaxies (as a mid-way between SO and Sc galaxies) con¬ 
structed by Bicker et al. 1200411 according to their chem¬ 
ical evolution models. In these spectra, effects of evolu¬ 
tion and redshifting had already been taken into account 
(therefore the correction actually includes also evolution 
correction). Our share was to calculate the relation be¬ 
tween the observed flux and the redshifted standard U, 
V or B filter flux for each set of cameras, filters and red¬ 
shifts. The throughput curves of the NICMOS filters and 
STIS clear imaging were taken from NICMOS and STIS 
Instrument Handbooks, respectively. Due to the available 
choice of filters, the final k-corrections for both rest-frame 
B-band and (U — V ) color remained rather modest, typi¬ 
cally around 0.4 — 0.7 magnitudes. 
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The final, k-corrected surface brightness profiles in 
rest-frame .B-color with estimated error bars are presented 
in the upper panels of Fig [3 In the lower panels, k- 
corrected rest-frame (U — V) profiles and corresponding 
error bars are given. To calculate the extent of the error 
bars, inaccuracy of our Hubble type classification and the 
population synthesis models were taken into account. The 
latter was difficult to estimate precisely because no actual 
uncertainties of chemical evolution models are available. 
Thus we made a rough estimate on the basis of deviations 
between synthetic spectra calculated by different authors. 
These uncertainties are amplified by different widths and 
shapes of the filter passbands. The uncertainties of the 
synthetic spectra at higher redshifts are larger than at 
lower redshifts, but the estimated errors remain compara¬ 
ble (0.25 — 0.45 mag). This is because of the spectra being 
stretched as (1 + z) with respect to the filter passband on 
one hand, thus enabling more exact determination of the 
flux, and the usage of the wide STIS imaging throughput 
at the tricky [/F-region for lower redshift galaxies on the 
other hand, which can not be very accurate. The step of 
determining actual flux according to the synthetic spectra 
dominates in the estimates of the uncertainties of color 
measurements on Fig Q 

The usage of different filters and detectors for imag¬ 
ing at different redshifts would cause errors in absolute 
photometry and introduce artificial trends with respect 
to redshift, if calibrations or k-corrections were calculated 
improperly. Considerable errors of this kind would be seen 
as a jump in magnitude-redshift plots at the redshift of the 
filter change, e.g. between z = 1.1 and z — 1.4 in our case 
of (U — V) vs. z graph in Fig0 The graph of absolute B 
magnitudes Mb vs. z was also checked for this effect; no 
jumps were found. 


4. Fitting by model profiles 

Although in some cases (usually at smaller redshifts), the 
bulge component was distinguishable from the disk, no 
attempt was made to split the luminosity profile into two 
components, because this would have given systematically 
different disk parameters at different redshifts. 

Luminosity profiles were fitted by a Sersic surface den¬ 
sity distribution (Sersic. 119681 ): 

1 (a) = /(0)exp[-6 n (r/r e ) 1/Tt ], (1) 

In this formula, r is the distance along the galactic ma¬ 
jor axis, r e is the effective radius, containing half of the 
total luminosity, n is the shape parameter, setting the 
curvature of the profile, and b n is a normalizing con¬ 
stant, dependent on n, enabling to keep r e the half-light 
radius. To determine 6 ra , either gamma functions or in¬ 
tegral equations have to be numerically solved. Usually, 
b n = 1.9992n — 0.3271 (Capaccioli 119891 ) or other formu¬ 
lae are used (e.g. Prugniel & Sirnien [l"997l Moriondo et 
al. flDiM Ciotti & Bert in HM5 l. which all can be success¬ 
fully applied at n values above n ss 0.5. If lower n values 



n 


Fig. 3. Normalizing parameter b n in Sersic formula CO 
at low n values. Comparison of our approximation o 
with numerically calculated actual values, together with 
approximations by Capaccioli 1111591 and Ciotti & Bertin 

nm 

are present, these solutions start to mislead. Therefore, we 
have calculated our own approximation, 

b n = 2n — — + (2) 

which can well be used down to n ~ 0.1. Comparison 
of our approximation 0 with numerically calculated val¬ 
ues together with approximations by Capaccioli 119591 and 
Ciotti & Bertin !T5i3§l is given in Fig. []] To emphasize de¬ 
viations at small n values, logarithmic scale is used. 

As a tool for fitting surface brightness distribution, the 
Sersic law has a purely empiric background. Starting from 
a spatial density distribution law and projecting it along 
the line of sight would be a more physical method, also al¬ 
lowing a straightforward comparison with kinematic data 
and thereby enabling the construction of self-consistent 
models for mass distribution. Thus, fitting of a surface 
density distribution, deduced from the space density dis¬ 
tribution law of the very general form 

p(a) = p( 0) exp [—(a/(fco 0 )) 1/n ], (3) 

was also tested for comparison. Here, p(0) is the central 
density, a = \JR 2 + z 2 /q 2 , where R and z are two cylin¬ 
drical coordinates and q is the axial ratio; ao is the char¬ 
acteristic radius and k is a normalizing parameter, similar 
to Sersic b n (see Einasto & Haud fT3S51 and Tenjes et al. 
119941 for more details about this distribution and reason¬ 
ing for its usage). This method gives the luminosity profile 
a somewhat different shape, but if fitted to a real density 
distribution, the main parameters (n, ao ss r e ) remain 
similar to those acquired from the Sersic law. 

Comparison of these two distributions at different n 
values is given in Fig 0] Surface brightness value p, and 
radius r are given in the same relative units for both dis¬ 
tributions. We stress that no rescaling has been conducted 
for matching the distributions; the fit has been achieved by 
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Fig. 4. Comparison of the Sersic surface brightness distri¬ 
bution (solid lines) and a distribution, derived from the 
space density distribution of Eq. 0 (dashed lines) at dif¬ 
ferent n values. Surface brightness value p and radius r 
are given in the same relative units for both distributions. 
No rescaling has been conducted to match the curves. 


keeping the total luminosities equal. Note that a perfectly 
exponential distribution, corresponding to the Sersic law 
with n=l cannot be attained with distribution 0 . 

For galaxies in the present sample, no additional kine¬ 
matic data are available at present. For this reason, we 
only present here the parameters resulting from fitting of 
the Sersic law 0 to the brightness profiles. This also al¬ 
lows more direct comparison with other studies. 


5. Results and discussion 

The rest-frame B luminosity profiles, fitted to the Sersic 
law, and the (U — V) color profiles with linear fits are 
shown in Fig|2] The main properties of the galaxies and 
the results of the photometric analysis are given in Tablc0 
Distribution of the parameter n, the effective radius r e , 
the (U — V) color distribution, and color gradient A (JJ — 
V)/Ar are shown in Figs. 0] [7] as a function of redshift. 
The three different brightness levels of the data points 
indicate three absolute rest-frame B magnitude-bins of the 
galaxies: white points stand for galaxies with Mb < —20 
mag, grey points for Mb = —19... — 20 mag and black 
points for Mb > —19 mag. For equal treatment of the 
galaxies at different redshifts, absolute brightness has been 
calculated within rest-frame surface brightness pb < 24 
mag/arcsec 2 . 



z 


Fig. 5. The Sersic luminosity profile curvature n as a func¬ 
tion of redshift. Open circles stand for galaxies with abso¬ 
lute luminosity Mb < —20 mag, gray circles for —19 mag 
< Mb < —20 mag and black circles for Mb < —19 mag. 



z 


Fig. 6. The Sersic luminosity profile effective radius r e 
as a function of redshift. Open circles stand for galaxies 
with absolute luminosity Mb < —20 mag, gray circles 
for —19 mag < Mb < —20 mag and black circles for 
Mb < —19 mag. The solid line shows a least-squares linear 
fit to the current sample r e (z)/r e ( 0) = 1 — 0.26z; the long- 
dashed line is the relation r(z)/r( 0) = 1 — 0.27z derived 
by Bouwens & Silk (120021) : the dotted and short-dashed 
lines give the relations r e ~ H(z)~ 1 and r e ~ H~ 2 ^ 3 (z), 
respectively. 

5.1. The shape of the luminosity profiles 

As is evident from the present sample, nearly half of galax¬ 
ies cannot be fit to an exponential disk model, while for the 
rest, n ~ 1 works well out to the highest redshifts (Fig.0. 
A similar result was obtained by Moth & Elston (120021) 
for the HDF North field. As suggested in Tamm & Tenjes 
(120051) . the more curved profiles might indicate the effect 
known as disk truncation, discovered by studies of local 
edge-on disk galaxies (de Grijs et al. 3()() 1 Pohlen et al. 
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Table 2. General galactic parameters. 


Name 11 

NED name 

(also J2000 coordinate) 


Mg 

(mag) 

Lb 

10 1O L(T) 

(kpc) 

n c 

U-V 

(mag) 

A(U-V)/Ar 

(mag/kpc) 

489 - 939 

HDFS 

J223252.31-603827.3 

0.49 

-16.4 

4.3 

1.8 

0.8 

1.55 

0.03 

607 - 767 

HDFS 

J223251.11-603840.2 

0.51 

-17.7 

3.3 

2.3 

1 

1.40 

-0.05 

381 - 481 

HDFS 

J223253.42-603901.6 

0.6 

-18.8 

3.2 

6.7 

0.7 

1.66 

-0.01 

708 - 695 

HDFS 

J223250.08-603845.6 

0.66 

-19.3 

3.2 

4.5 

1 

1.26 

0.02 

317 - 817 

HDFS 

J223253.65-603846.9 

0.89 

-18.3 

1.9 

3.0 

1.2 

1.28 

-0.01 

380 - 1027 

HDFS 

J223253.43-603820.6 

0.99 

-20.8 

1.8 

4.1 

0.7 

- 

- 

419 - 572 

HDFS 

J223253.03-603854.8 

1.01 

-20.1 

1.7 

4.8 

0.75 

1.69 

-0.10 

648 - 814 

HDFS 

J223250.70-603836.6 

1.12 

-18.6 

1.6 

1.3 

0.8 

1.66 

0.08 

836 - 1023 

HDFS 

J223248.78-603821.0 

1.37 

-21.1 

1.3 

3.5 

0.75 

1.18 

-0.01 

854 - 606 

HDFS 

J223248.59-603852.2 

1.37 

-19.6 

1.2 

0.7 

0.6 

1.05 

-0.02 

353 - 444 

HDFS 

J223253.71-603904.4 

1.39 

-20.1 

1.2 

1.9 

1 

1.22 

0.02 

507 - 451 

HDFS 

J223252.13-603903.9 

1.65 

-20.2 

1.0 

2.0 

1.1 

1.14 

-0.01 

786 - 619 

HDFS 

J223249.29-603851.2 

1.65 

-17.2 

0.8 

1.4 

0.3 

1.41 

0.00 

472 - 850 

HDFS 

J223252.49-603833.9 

1.89 

-19.3 

0.8 

2.0 

1.1 

0.95 

0.09 

596 - 1002 

HDFS 

J223251.22-603822.6 

1.93 

-20.0 

0.8 

2.6 

0.6 

1.23 

0.01 

647 - 899 

HDFS 

J223250.70-603830.3 

2.01 

-19.7 

0.8 

2.8 

1.1 

1.71 

-0.06 

565 - 941 

HDFS 

J223251.54-603827.1 

2.33 

-19.7 

0.5 

1.6 

0.6 

0.64 

0.25 

501 - 899 

HDFS 

J223252.19-603830.2 

2.39 

-20.8 

0.4 

1.7 

0.9 

0.64 

0.05 

891 - 961 

HDFS 

J223248.22-603825.6 

2.39 

-20.8 

0.3 

1.0 

1 

0.61 

0.00 

899 - 591 

HDFS 

J223248.13-603853.4 

2.55 

-19.3 

0.2 

0.9 

0.9 

0.62 

0.00 

784 - 511 

HDFS 

J223249.30-603859.4 

2.56 

-19.2 

0.1 

1.5 

1 

0.58 

0.05 

651 - 392 

HDFS 

J223250.67-603908.3 

2.58 

-19.8 

0.1 

2.7 

0.5 

0.59 

0.20 


“ The x and y pixel coordinate on HDF-S NICMOS final images 
b Photometric redshifts are from Yahata et al. BSp. 
c Steepness parameter of the best-fit Sersic profile. 


6002 ) and more recently detected also at higher redshifts 
(Perez f?0l) 1 Trujillo & Pohlen 2()()■' )•. Such a luminosity 
distribution is sometimes fitted by a double-exponential 
profile (Pohlen et al. BD'iM Perez) 2(j(T4 ). The physical back¬ 
ground of disk truncation has not become clear as yet (see 
e.g. Sasaki ITM7I and de Grijs et al. 2001 ] for more discus¬ 
sion about the origin of this phenomenon). 

In contrast to these truncated disks, some galaxies in 
the present sample have flattening profiles (e.g. 607 — 767, 
708 — 695, 380 — 1027, 353—444). Similar profiles have been 
measured near z =1 for gannna-ray burst-selected disks 
(Conselice et al. 12005bi and out to z =2> for large disks 
in HDF-South primary field (Labbe et al. 1200311 . Erwin 
et al. m have found a significant proportion of lo¬ 
cal barred SO-Sb galaxies with flattening profiles, which 
they call “anti-truncated” disks. They claim that, at least 
in several cases, this effect might be caused by interac¬ 
tions. Observations of interacting galaxies (Chitre & Jog 
12002 1 and numeric simulations (Bournaud et al. 2ti0'5 l 
have demonstrated the possibility of mergers to smear out 
the outer parts of luminosity profiles. In the case of our 
sample, visual companions are present for all four galax¬ 
ies with flattened profiles, thus interactions may have a 
role in the development of this flattening. The compan¬ 
ions have been masked prior to ellipse fitting, therefore 
the luminosity of the companions cannot be responsible 
for the effect. 


The 7 galaxies studied in Tamm & Tenjes 125031 20()5 ) 
have significantly low n values compared to the present 
sample. This is most likely a result of the small size of 
the samples or selection effects, favoring luminous galaxies 
with regular kinematics. 

5.2. Sizes of disks r e 

In Fig. EO evolution of the size of the present sample is 
shown. Error bars of sizes were derived by varying r e while 
fitting the Sersic law within the surface brightness profile 
error bars. The slight degeneracy between parameters n 
and r e has also been taken into account. 

Starting from the locally observed Universe and mod¬ 
eling galactic evolution backwards, Bouwens & Silk ■ 
derived a scaling relation r( 0 )(r(O) = 1 — 0.270 for the 13- 
band radius. The long-dashed line in Fig. EJ shows this 
relation. The least-squares linear fit to our sample gives 
r e {z)/r e { 0) = 1 — 0.260 (solid line). Here, the median value 
of Sersic half-light radii ( r e ) ss 4.2 kpc for late-type galax¬ 
ies from the SDSS survey (Shen et al. 120031) has been used 
as the value of r e (0) (considering the median absolute lu¬ 
minosity of the present sample). It is seen that the model 
developed by Bouwens & Silk ■ fits our data rather 
closely. Our result is consistent with the small r e values 
found for UDF galaxies by Elmegreen et al. (120051) . 
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Fig. 7. Galactic parameters - rest frame (U — V ) color 
(a) and A (U — V)/Ar (b) - as a function of redshift. 
Designations are the same as in Fig. 0 



Fig. 8. Effective surface brightness g e , plotted as a func¬ 
tion of log[r e ). Designations are the same as in Figs. 00 
Solid line shows a least-square fit. 

Galactic size evolution has also been predicted from 
simple models of the hierarchical structure formation sce¬ 
nario and an approximate scaling has been derived with 
disk formation time: r e ~ V v i T /H(z), where in the case of 
D 0 = 1, H{z) = + ~) 3 + Da] 1/2 (Mo et al. 11998)) . 


In a study of galaxies at a redshift range z = 1.5 — 5, 
Ferguson et al. mm found the scaling of r e ~ H 1 (z) 
to give a good fit to their sample, while Trujillo et al. 
m achieved a better fit with r e ~ H 2 / 3 (z) for galax¬ 
ies at z = 0.3 — 2.5. The relations r e ~ H~ 1 {z) and 
r e ~ H~ 2 / 3 (z), scaled with r e (0) « 4 kpc, are also shown 
in Fig. 0 as dotted and short-dashed lines, respectively. 
The former shows a decrease significantly faster than that 
of the present sample. However, it must be kept in mind 
that the theoretical relation uses the disk formation red¬ 
shift instead of the observed redshift for scaling. Further 
more, in contrast to the present sample and the literature 
referred to above, Ravindranath et al mm and Cassata 
et al (120051) did not detect any significant evolution of 
galaxy sizes with redshift in a thorough study of galaxies 
at z < 1.2. 

Fig. 0 shows that the trendline fitted to the present 
sample is quite sensitive to the data point corresponding 
to the largest galaxy with r e = 6.7 kpc at z = 0.6. Its 
extensive disk is obvious also on Fig.0and its light profile 
is a very regular one, thus there is no doubt that this 
galaxy is clearly a large disk galaxy. The only reason for 
this data point to be significantly offset could be a wrong 
photometric redshift estimate. 

5.3. Color (U — V) and color gradient A(U — V)/Ar 

In Fig. 0U rest frame (U — V) color as a function of red¬ 
shift is presented. The (U — V) color shows a mild decrease 
until z ~ 2, followed by a notable drop of roughly 0.5 mag 
by 2 « 2.5. At lower redshifts, the dependence can be 
compared to the Deep Groth sample studied by Weiner 
at al. (120051) . for which a very slight, if any, decrease of 
(U — B) was detected for late-type galaxies in the redshift 
range 2 = 0 — 1.5. The notable drop of the (U — V) val¬ 
ues at 2 « 2 cannot be due to miscalibration - a uniform 
set of filters and calibrations has been used for the red¬ 
shift range 2 = 1.4 — 2.6 (see Sects 2 and 5.5). The drop 
can rather be related to a major star formation peak at 
2 > 2 (see below). A decrease of the rest-frame (U — V) 
color with increasing redshift up to 2 ~ 3 was derived 
also by Kajisawa & Yamada mm . More detailed com¬ 
parison with our results is not possible at present because 
in their study, different morphological types are presented 
together. Our result is quantitatively close also to the re¬ 
sult obtained from the analysis of the synthetic spectra of 
Bicker et al. mm . Absolute (U — V) values are rather 
sensitive to k-correction, which, in turn, is uncertain. This 
may cause a constant vertical shift of (U — V) values. 
However, trends in Fig. 0 should be independent of these 
uncertainties. 

Fig- 0^> shows that at redshifts 2 < 2, there are small 
or no color gradients. Gradients begin at 2 > 2 and are 
dominantly positive. A similar rise with redshift of the 
rest-frame color (C /218 — U 300 ) gradients was discovered by 
Moth & Elston 1)20021 . suggesting that star-formation was 
more centrally concentrated at higher redshifts. 
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Recent hydrodynamical simulations of galaxies in 
ACDM cosmology by Robertson et al. mm indicate that 
star formation in disks peaks between redshifts z = 2 — 4. 
Studies of the Fundamental Plane of early-type galaxies 
at redshifts 2 < 1.3 have also shown that the last epoch 
of major star formation peaks at z ~ 2 — 3.5 (see Holden 
et al. 120051 and references therein). A jump in star forma¬ 
tion beyond redshift z = 2 matches well with the jump in 
(U — V) color in Fig.Qi (see above). 

Smaller color gradients at lower redshifts z = 0.5 — 2 
are due to smoothing by interactions in disk evolution 
(Conselice et al. l2005aJI and large scale gas motions 
(Samland & Gerhard. 120031 ) . 

Positive color gradients also exist in local disk galaxies 
(MacArthur et al. ; 2'(')T) 1). Positive color gradients may in¬ 
dicate disk formation from inside out. For the Milky Way 
protogalaxy such formation scenario was suggested by van 
den Bergh mm- In the case of high-redshift galaxies, di¬ 
rect interpretation of color gradients may be rather com¬ 
plicated. Color gradients are determined by radial distri¬ 
bution of initial metallicities, stellar ages, star-formation 
rate, gas accretion details etc. Dependently on model de¬ 
tails, both outside-in and inside-out formation of disks can 
be simulated (Sommer-Larsen et nl. 2003 Robertson et al. 

mm- 

According to the analysis of their disk galaxy forma¬ 
tion model, Westera et al. m have found that lumi¬ 
nosity from U to V mainly indicates star formation, and 
the effects of age and metallicity are negligible. Thus the 
positive color gradient most possibly refers to intensive 
star formation in central parts of galaxies at z > 2. 

Relying on (U — V), we may roughly estimate the disk 
formation time. According to our photometry, the disks 
at z ~ 2.5 have (U — V) ~ 0.7. On the basis of simple stel¬ 
lar population chemical evolution models (Worthev il9941 
Bressan et al. mg, these disks have ages 1 — 2 Gyr and 
have thus formed at z = 3.5 — 7. 

5.4. Surface brightness at effective radius 

Coenda et al. (120051) have shown that disks obey similar 
photometric scaling relations as spheroids and elliptical 
galaxies. Degeneracy exists between the effective radius, 
effective surface brightness, the parameter n and abso¬ 
lute magnitude. For the current sample, a trend can be 
shown for the effective surface brightness /i e to diminish 
with larger effective radius r e , presented in Fig. 0 Least 
squares fit gives g e = 3.0 log r e + 19.4 mag, which, re¬ 
garding the small size of the sample, is remarkably close 
to the dependence found by Coenda et al. (120051) g e = 
3.4 log r e + 18.7 mag . 

5.5. Uncertainties of the study 

Small sample sizes and selectional biases are critical top¬ 
ics for most of the studies of high-redshift galaxies. The 


present sample of 22 galaxies can not offer reliable statis¬ 
tics and cautions against jumping into severe conclusions. 

Morphological identification of these disk galaxies was 
obtained by visual inspection of 77-band exposures, all 
galaxies suspected for being disks were included in the 
sample. In some cases, also the high-resolution STIS 14- 
band images were checked. Some of the galaxies of the 
present sample have been studied on the basis of NICMOS 
images by Rodighiero et al. and classified as E- 

S0. We would like to point out that for classification and 
accurate photometry it is important to carefully decon¬ 
volve the PSF of NICMOS images. Even in the case of 
good restoration of NICMOS images, additional observa¬ 
tions with other detectors, preferably with higher resolu¬ 
tion (e.g. STIS in the present case) still prove beneficial. 
For example, visual inspection of the STIS image of the 
galaxy 419 — 572 reveals its nature of a late-type (at least 
Sb) galaxy. 

Usually, the cosmological dimming as (1 + z) A favors 
the selection of more luminous galaxies at higher redshifts. 
For our sample, Figs 0 El reveal no clear dependence of 
any of the parameters on the absolute luminosity of the 
galaxies; this fact allows us to consider the sample to be 
sufficiently complete for the conducted analysis. 

The selection of the sample on the basis of 77-band 
imaging (corresponding to 77-band at 2 > 2) might cause 
a slight bias towards bluer galaxies at higher redshifts. On 
the other hand, the change of the spectral energy distribu¬ 
tion is slow between B and R passbands and we consider 
this possible selection effect not to be a reason of a color 
jump at 2 ~ 2 in Fig.0 

For calculating k-corrections, all galaxies were as¬ 
sumed to be Sb morphological types. This approximation 
causes additional scatter of absolute magnitudes and col¬ 
ors. However, according to our estimates, the resulting 
uncertainties remain within 0.2 magnitudes even in the 
worst cases. 

Among possible sources of errors and uncertainties, 
the ever insufficient imaging depth and spatial resolution 
should be considered. With the currently available NIR 
imaging a dilemma remains, whether to consider faint 
patches around galaxies to belong to the galaxy and in¬ 
clude them in photometry or to mask them as extrane¬ 
ous. This may slightly affect the luminosity distribution 
at the outer regions of a given galaxy and the values of 
the photometric parameters, especially the shape param¬ 
eter n . In the present case, we have followed a rather 
conservative masking strategy, excluding most of the “sus¬ 
picious” patches. A more liberal treatment would give 
slightly higher n and r e values and lower absolute magni¬ 
tudes. The shift of these values would be rather system¬ 
atic. 

Insufficient spatial resolution influences the interpre¬ 
tation of surface brightness profiles in the central parts of 
galaxies. In the present study we decided not to decompose 
galaxies into the bulge and disk components. Otherwise, 
rather serious systematic differences would appear in han¬ 
dling of low-redshift and high-redshift galaxies. In general, 
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decomposing light profiles into a share of a bulge and a 
disk would decrease the disk radius. 
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